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ABSTRACT 



It is generally assumed that the magnetic fields of millisecond pulsars (MSPs) 
are ~ 10 8 G. We argue that this may not be true and the fields may be appreciably 
greater. We present six evidences for this: (1) The ~ 10 s G field estimate is based 
on magnetic dipole emission losses which is shown to be questionable; (2) The 
MSPs in low mass X-ray binaries (LMXBs) are claimed to have < 10 11 G on 
the basis of a Rayleygh- Taylor instability accretion argument. We show that 
the accretion argument is questionable and the upper limit 10 11 G may be much 
higher; (3) Low magnetic field neutron stars have difficulty being produced in 
LMXBs; (4) MSPs may still be accreting indicating a much higher magnetic field; 
(5) The data that predict ~ 10 8 G for MSPs also predict ages on the order of, and 
greater than, ten billion years, which is much greater than normal pulsars. If the 
predicted ages are wrong, most likely the predicted ~ 10 8 G fields of MSPs are 
wrong; (6) When magnetic fields are measured directly with cyclotron lines in X- 
ray binaries, fields 3> 10 8 G are indicated. Other scenarios should be investigated. 
One such scenario is the following. Over 85% of MSPs are confirmed members 
of a binary. It is possible that all MSPs are in large separation binaries having 
magnetic fields > 10 8 G with their magnetic dipole emission being balanced by 
low level accretion from their companions. 

Subject headings: magnetic fields, millisecond pulsars 
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Introduction 



The discovery of millisecond pulsars (MSPs) was reported by IB acker et al.l ( 119821 ). 
'heir origin was su ggested to be due to spin-up by accretion torques of a binary companion 



( lAlpar et al. 



19821 ) such as in low mass X-ray binaries (LMXB). The strength of the 
magnetic field of the MSPs is generally estimated assuming that they are isolated and there 
is no accretion. Only energy loss due to magnetic dipole emission is assumed. From the 
spin-down rate, the magnetic field is estimated to be ~ 10 8 G, which is universally used. All 
of the articles published on the MSP since 1982 have assumed the validity of this scenario 
and the relatively low magnetic field ~ 10 8 G. This can be compared with the typical 
~ 10 12 G field of pulsars. In this article we argue that the magnetic fields of MSPs may be 
> 10 8 G. 

Assuming that the spin-down of a pulsar is due to magnetic dipole emission, the surface 
dipole magnetic field is obtained from the relation 

B = 3.2x 10 19 (PP) 1/2 G, (1) 



where the neutron star moment of inertia is assumed to be 10 g cm , t 



re neutron star 



radiu s 10 16 cm, P is the pulsar period in seconds, and P is its derivative ([Manchester et al. 



2005|). 



MSPs are characterized by spin-down rates four to six orders of magnitude less than 
normal pulsars implying ages ~ 10 10 years. It is argued that they probably require their 
short periods through a recycling process in which mass and angular mom entum are 



trans ferred to an old and slowly rotating pulsar from a binary companion ([Manchester 
200J). 



We discuss in section II "Can the magnetic field of MSPs be determined assuming 
magnetic dipole emission?" , in section III "Can the magnetic field of a MSP be determined 
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by the rate of accretion in LMXBs?", in section IV "Can low magnetic field MSPs be 
produced in LMXBs?", in section V "Are MSPs still accreting?", in section VI "Are the 
indicated ages of MSPs realistic?" and in section VII "Measurement of magnetic fields in 
X-ray binaries with cyclotron lines". Our "Conclusions and Discussion" are presented in 
section VIII. 



2. Can the magnetic field of MSPs be determined assuming magnetic dipole 



emission : 



In most theoretical models of pulsars the braking torque is proportional to some power, 
n, of the rotation frequency Q (= 2n / P) 



n 



-KQ n 



(2) 



where Q is the time derivative of Q ( [Manchester fc Taylorl 119771 ). The parameter n is 
known as the braking index and K is a positive constant. For braking by magnetic dipole 
radiation or particle acceler ation in a dipole field, which is generally assumed, n = 3 



( [Manchester Sz Taylor 



19771 ) 



The very low spin-down rates of the MSPs have so far precluded any direct 
measurements of the braking index. An index of n = 3 is indicat ed for old radio pulsars, 



but values n ~ 1.5 — 2.8 have been mea sured in younger pulsars (ILyne 



2004 



Ferrario fc Wickramasinghe 



20070 . 



1996 



Hobbs et al. 



Pacini 



( 119671 ) pointed out that a rotating magnetized star would radiate, in the form of 
electromagnetic waves, magnetic dipole radiation at the rotation frequency. The radiation 
reaction torque transmitted to the star by the magnetic field is 

2(/i sin a" 2 



N 



3c 3 



(3) 
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where \x is th e magnetic-clippie mome nt an d a is the angle betwee n the magnetic and 



rotation axis. 



Gunn fc Ostrikerl (119691 ) and 



Ostriker fc Gunnl (119691 ) recognized that this 



torque could account for the observed secular increase in pulsar periods. Since N oc Q 3 , the 
braking index n is three for these pulsars. Since \x ~ BqR 3 , where Bq is the surface magnetic 
field and R is the radius of the star, the surface magnetic field Bq can be estimated from 
the observed period derivative (with sin a assumed approximately equal to one) 



Bn 



11/2 



(4) 



where I is the moment of inertia. For I ^ 10 45 gem 2 , R = 10 e cm and P in seconds we 



obtain Eq. (1) ( ([Manchester fc Taylor 



1913)) 



It is not clear how an isolated MSP is formed, but it is reasonable to assume that it is 
not found in a vacuum. A low density plasma undoubtedly exists around the MSP. In that 
case, the electromagnetic dispersion relation has the form 



2 ;„2 



[e.g., 



Sturrock 



19941 ) where v p , e (= u P)e /27i) is the plasma frequency, 



(5) 



(6) 



and n e (in units of electrons per cm 3 ) is the electron density. For u < u pe , the wavenumber 
k is imaginary and the electromagnetic wave is rapidly damped. Thus an electromagnetic 
wave of frequency v = uj/2h = 10 3 , emitted by a MSP, will not propagate away from it if 
the density of the plasma is n e > 0.013. This limit for the plasma density is very low. It can 
be compared to n ~ 1 for the interstellar medium in our galaxy and n > 1 in a globular 
cluster. 

Particle acceleration by the rotating magnetic field of the neutron star was studied by 



Goldreich &: Julian! ( 119691 ). The rotation of the neutron star and its associated magnetic 
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field produces strong electric fields in the space surrounding the star. iGoldreich Sz Julian 
were the first to point out that, because of these fields, the region surrounding the star 



cann ot be a vacuum and must contain a substantial space charge ([Manchester Sz Taylor 



1973). 



The work of IGoldreich Sz Julian! however, considered a simplified system (an 
axisymmetric rotator) whereas in fact pulsars must possess a non-axisymmetric field (an 
oblique rotator) in order to produce a periodic signal. Furthermore, they did not give a 
self-consistent description of the currents and fields surrounding the star. Their model 
involves flow of charge of one sign through a regi on of space with charge of the opposite 



sign, a situation unlikely to exist in a real pulsar ( Manchester Sz Taylor 



1973). 



From the above we can not be certain that the MSPs are radiating as a magnetic 
dipole and that the P and P of MSPs give a true value of B . 



3. Can the magnetic field of a MSP be determined by the rate of accretion in 

LMXBs? 



It has also been argued that MSPs in LMXBs have low magnetic fields, < 10 11 
G, due to the accretion pressure produced by a weak Ray leigh- Taylor instability. The 



argument is that, if the field is hi gher, accretion onto t 



re neutron star could not occur 



and produce the observed X-ravs (Spruit Sz Taam 



Psaltis Sz Chakrabarty 



1999 



Lamb 



1990 



Miller. Lamb, Sz Psaltis 



1998 



20061 ). If the stellar magnetic field is > 10 n G, all of 



the gas in the disk will couple to the neutron star and be funneled out of the disk plane, 
toward the magnetic poles of the star. For magnetic fields < 10 11 G, some of the accreting 
gas will continue as a disk flow towards the stellar surface as a result of the Rayleigh- Taylor 
instability and produce the observed X-rays. However, if the instability is greater than the 
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Rayleigh- Taylor instability, magnetic fields > 10 G are possible. 

Previous analysis neglected the strongest disk instability, which is m uch greater than 



the Rayleig h- Taylor instability: the magneto-rotational instability (MRI) (IBalbus &: Hawley 



1991 



19981 ) . In a Keplerian disk the maximum growth rate of the MRI is ~ (3/4)0, with a 



typical perturbation wavenumber Q/Va, where O(Va) is the angular (Alfven) velocity. The 



disk becomes very dynamic and ful l of reconnection events f lKudoh. Matsumoto. fc Shibata 



2002 



Machida fc Matsumoto 



20031 ). 



We thus conclude that the low accretion rate in LMXB could occur for a stellar 
magnetic field > 10 11 G due to the MRI. 

We discussed above one of the most popular scenarios for determining the magnetic 
field of a MSP by the rate of accretion in LMXBs. Other means of accretion are possible. 
In the magnetic Cataclysmic Variables, where a magnetic white dwarf accretes mass from a 
low-mass companion star, the existence and extent of the disk strongly correlate with the 
strength of the magnetic field of the white dwarf. Hence, in the most strongly magnetic 
system (10 7 — 10 8 Gauss: the AM Herculis systems), accretion occurs directly via a funnel, 
without the formation of a disk. These strong fields are measured directly via the detection 
of strong cyclotron lines and/or Zeeman lines from the surface of the white dwarfs. In the 
Intermediate Polars, accretion occurs via a truncated disk, the extent of which is related 
to the strength of the field of the white dwarf. In some cases, fields in the Intermediate 
polars have also been measured, directly, via the detection of polarized radiation and, more 
recently, of cyclotron lines in the polarization spectra. 
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4. Can the low magnetic field MSPs be produced in LMXBs? 



It has been suggested that pulsa r magnetic fields are lowered by recycling in binaries 



see 



Bhattacharya fc van den Heuvell 1991 for a review). The mechanism for this remains 



uncle ar, with suggestions including decay of crustal fields due t o heating ( 



19861 ). burial of the field (jBisnovatyi-Kogan fc Komberg 



1974 



Romani 



Blondin fc Freese 



1990 



decay o 



19901)). 



core fields due to flux tube ex pulsion from the superfluid interior ( jSrinivasan et al. 



19931) and 



Konar &: Bhattacharya! (119971 ) suggested that rapid ohmic decay in the accretion 



heated crust occurs. On the one hand, the heating reduces the electrical conductivity and 
consequently the ohmic decay time-scale induces a faster decay of the field. On the other 
hand, the deposition of matter on top of the crust pushes the original current carrying 
l ayers into deeper and denser regions where the higher conductivity slows down the decay 



( jKonar &: Bhattacharya 



1993). 



In a class of recycling models in which the magnetic field decrease is a function only 
of the amount accreted onto the neutron star, it was shown that no model of this class is 



consistent with all available data (IWijers 



19970 



The detection of coher ent X-ray pulsations with a millisecond period in a handful of 

ten u sed in support of the idea of accretion-induced field 



LMXB s (lLamb & Yu 



2005 ) is o 



decay (jWijnands &: van der Klis 



19981 1 . However, whether this is evidence simply for field 



submersion and spin-up d uring the accretion disk phase, or for field decay and spin-up 



remains to be established (jFerrario Sz Wickramasinghe 



20071 ) 



We conclude from the above, that it is not clear if low field ~ 10 8 G field MSPs can be 
produced in LMXBs. 

We pointed out that the field could simply be s ubmerged by the ac creted matter and 



would then re-emerge later on when accretion stops. 



Zhang et al. 



(120091 ) show that there is 
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no evidence for field restructuring and/or decay in accreting magnetic white dwarfs. This 
fact was already known for the isolated magnetic white dwarfs. As far as we know, this 
could apply also to neutron stars. 



5. Are MSPs still accreting? 

Almost all (> 85%) MSPs have been confirmed to be in binaries. The central idea 
behind the MSP-LMXB connection is that LMXBs can provide the long-lived phase of 
moderate mass transfer rates thought necessary to spin-up the neutron star to millisecond 
periods. The MSP is supposed to be produced when the accretion shuts off. 

But does the accretion ever shut-off? In order to have a shut-off of the mass accretion 
rate, M, the accretion evolution time 

m = w (?) 

needs to be shorter than the spin-down time scale 

r spm = oc M- 3 / 7 (8) 

*spin 

or 

^ 7~spin- (9) 

The accretion rate, M, decreases as the orbital separation grows. At late times, for 
gravitational wave driven orbital evolution in the limit when the low mass companion of 
the neutron star is very much less than that of the neutron star, 

m oc M- 11 / 14 (10) 
Thus as M decreases, the ratio T sp i n /T^ decreases. For a neutron star with a magnetic field 
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B = 5 x 10 8 G, the two timescales are roughly equal at M ~ 10 9 Mp)/yr. Th e system thus 



has time to come into equilibrium and accretion never shuts off (IDeloye 



20081 ). 



It was noted by 



Deloyd ( 120081 ) that in two of the four classes of LMXBs, mass transfer 



to the neutron star never shuts off. Both of these systems ha ve gradu ally declining mass 



transfer rates, as the system evolves to longer orbital periods. 
MSPs can be associated with one of these two classes. 



Deloyd noted that all known 



Thus, MSPs may always be subject to a very low accretion rate. 



6. Are the indicated ages of MSPs realistic? 



In normal pulsars P and P are distributed about mean values P ~ 0.1s and 
P ~ 10- B ss -1 . This implies fields B ~ 10 11 - 10 13 G. Values of MSPs of P ~ 5ms and 
P ~ lO-^ss -1 imply fields of B ~ 10 s - 10 9 G. The age of the MSP is obtained from the 
relation 

P / P \ n ~ 1 " 

(11) 



T 



(n - l)P 

where Pq is the initial period. For P -C P and n = 3 



T 



P 

2P 



(12) 



Fr om this expression, however, the ages of many MSPs are more than ten billion 



years ( 



Ferrario fc Wickramasinghe 



bscano et al. 



20071 ). Half of the corrected sample of MSPs of 



(119991) have charact eristic ages comparable or greater than ten billion years 



(IFerrario fc Wickramasinghe 



20071 ). 



These ages are orders of magnitude greater than the ages of normal pulsars. Since the 
ages of the MSPs are obtained from P and P assuming the scenario that magnetic dipole 
emission is the only energy loss mechanism, and there is no accretion, a doubtful age makes 
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the scenario doubtful and the evaluation of Bq from P and P doubtful. 

7. Measurement of magnetic fields in X-ray binaries with cyclotron lines 

The only sure way of measuring magnetic fields is by observing the electron-cyclotron 
line. The many observed electron-cyclotron me asurements of magn et ic fields in X-ray 



binari es indicate fi elds > 10 8 G: 4U 03 52+309 



( 2001 ). GX 301-2 



Coburn et al. 



3) 



Coburn et al 



(120021), BA' 



fl2002h . B ATSE Pulsar t eam, 4U 1907+09 



SE Pulsar tea 



Coburn et al. 



3L 



(|2Q02|); 



Delgado-Marti et al. 



4U 1538-52 



Vela X-l Coburn et al. 



(120001 ). XTE 



4U 1626-67 



J19646+274 



Coburn et al 



Cen X-3 ICoburn et al 



Shira 



(I2002h 



( 200 



2): 



Coburn et al 



Bavkal et al 



(I2002h. BATSE P u 



Coburn et al 



2002 



sar team; A 0535+26 



(20 021), LMC X-4 



(200 



Negueruela et al. 



La Barb era et al 



van Kerkwiik et al 



(I2002h : iNelson et al 



kawa & Lai (2002); 4U 0115+63 



(1200 ll ) 



(119981 ) . BATSE Pulsar team, 



1997h. BATSE Pulsar tea m; V 0332+5 3 



Coburn et al 



Oosterbroek et al. 



([2002), Her X-l 



Coburn et al. 



20011 ). BATSE Pulsar team. No cyclotron measurement has 



indicated a magnetic field ~ 10 8 G on the neutron star. 

To be complete, the lack of detection of cyclotron lines in the X-ray spectra of LMXBs 
could be used, however, to support the view that their fields are low. The fields measured in 
accreting neutron stars using cyclotron lines are those of neutron stars in High Mass X-Ray 
Binaries, which are still very young and mostly wind accretors. If the fields of neutron stars 
in LMXRBs were of the order of 10 8 Gauss, the cyclotron lines would be in the optical 
spectral band. Such lines, arising from extremely small accretion shocks, would be swamped 
by the radiation from the accretion disk and would thus be hidden from detection. 



1 http:/ /www. batse.msfc.nasa.gov/batse/pulsar/index.html 



-12- 



8. Conclusions and Discussion 

The period P and its time derivative P are used to obtain the magnetic field Bq ~ 10 8 
G of MSPs. We presented six evidences in sections II- VII which indicate that the measured 
P and P may not indicate the true value of B and that B may be > 10 8 G. The evidences 
are: 

1. Magnetic dipole emission is assumed with a braking index n = 3. The braking 
index of MSPs, however, have never been measured. We show that the magnetic 
dipole electromagnetic emission does not propagate in realistic interstellar media and 
a self-consistent model of magnetic dipole particle emission has never been made, 
putting doubt on the necessity of the magnetic dipole emission assumption for MSPs; 

2. It is assumed that in LMXBs only the Rayleigh- Taylor instability occurs in the 
accretion disks and a magnetic field < 10 11 G needs to be present at the surface of the 
neutron star to allow for accretion to occur. The much stronger magneto-rotational 
instability is probably present allowing for a field ^> 10 11 G on the neutron stars; 

3. It has been assumed that the abnormal low magnetic fields ~ 10 s G for pulsars could 
be produced in LMXBs. There is, however, no self-consistent theory that does this; 

4. The standard scenario is that the pulsar is spun-up by accretion and then the 
accretion stops leaving the MSP to spin-down only by magnetic dipole emission. The 
accretion may never stop, however. A balance between dipole emission balanced by 
accretion could occur; 

5. The P and P of MSPs indicate ages > 10 10 years. This is comparable to the age of 
the Universe and very much older than other pulsars; and 

6. Measurements of cyclotron lines in X-ray binaries all indicate magnetic fields 3> 10 8 G. 
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Together, the six above evidences indicate that the general scenario that is used to 
deduce from P and P of MSPs that B ~ 10 8 G may be incorrect. We suggest investigating 
alternate scenarios. A possibility is that MSPs may have B 3> 10 8 G and that the observed 
P of MSPs is due to a balance of magnetic dipole emission with a low-level of accretion 
from a companion in a large separation binary. 

To be complete, the magnetic fields of MSPs, however, may be low due to the following 
scenario. The neutron stars were weakly magnetic at birth (maybe as a consequence of 
AIC). In this case, field decay would not be necessary. After all, only about 15% of white 
dwarfs are magnetic, with magnetic fluxes similar to those inferred for radio pulsars. It is 
not clear why one should assume that ALL neutron stars are highly magnetic. 

R.S.S. thanks the Brazilian agency FAPESP for financial support (2009/06770-2). 
R.O. thanks FAPESP (06/56213-9) and the Brazilian agency CNPq (300414/82-0) for 
partial support. We would like to thank an unknown referee for very helpful comments. His 
comments were included at the end of sections 3, 4, 6 and 8. 
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